2К 


УДК 658.325 


Д.С. Кныш, В.М. Курейчик 
ТИ ЮФУ в г. Таганроге, Российская Федерация 
ул5етап3 3 @уапдех.та, Каг@виге.га 


Параллельный генетический алгоритм 
с нечетким оператором миграции 


Развитие вычислительной техники способствовало развитию методов параллельных вычислений при 
решении оптимизационных задач. Генетические алгоритмы — эффективное средство решения различных 
оптимизационных задач. В этой связи необходимы стратегии применения параллельных вычислений 
в генетических алгоритмах. По аналогии с эволюционными процессами за основу параллельных 
вычислений в генетических алгоритмах взят механизм миграции. Авторами предлагается реализация 
оператора миграции с использованием аппарата нечеткой логики. Проведенные испытания разработанной 
модели показали ее эффективность. 


Введение 


В настоящее время вычислительные возможности компьютеров возросли за 
счет использования многоядерных процессоров. В связи с этим необходимо активно 
использовать механизмы распараллеливания для эффективного освоения потенциала 
вычислительных средств. Разработкой методов параллельных вычислений занима- 
ется междисциплинарное направление в науке — распределенные вычисления. В этой 
области проводятся исследования эффективности параллельных вычислений, фор- 
мируются стратегии распараллеливания и разрабатываются методы их реализации. 
Из всего разнообразия алгоритмов решения оптимизационных задач генетические 
алгоритмы, по аналогии с эволюционными процессами в природе, должны использо- 
вать параллельные вычисления для эффективного решения задач оптимизации [1-3], 
[4]. На протяжении вот уже двадцати лет продолжаются исследования использова- 
ния параллельных вычислений в генетических алгоритмах [5]. 

Генетические алгоритмы (ГА) базируются на модели биологической эволюции 
и методах случайного поиска [2], [3], [5-10]. ГА работают с популяцией индивидов 
(особей или хромосом), или альтернативных решений задачи [10]. Популяция разви- 
вается за счет применения механизма рекомбинации (объединение двух и более ро- 
дителей для получения потомков) и мутации (случайное изменение индивида). ГА 
эффективно применяются для решения задач оптимизации. Оптимизационную зада- 
чу можно описать следующим кортежем: <Х, О, О>, где Х является множеством 
всех решений (индивидов) для данной задачи, О — ограничения, накладываемые на 
множество Х, для получения допустимых решений и О — целевая функция, с по- 
мощью которой можно определить наилучшее (оптимальное) решение. 

Целью данной работы является разработка нового оператора миграции на ос- 
нове аппарата нечеткой логики для параллельного асинхронного генетического алго- 
ритма оптимизации многокритериальной функции. 

Статья построена следующим образом. В разделе 1 приведены основные меха- 
низмы распределения вычислений в ГА и их краткое описание. В разделе 2 представ- 
лена модифицированная модель параллельного генетического алгоритма. В разделе 3 
описывается нечеткий оператор миграции. В разделе 4 проведены исследования эф- 
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фективности работы модифицированного параллельного генетического алгоритма, ко- 
торые показали эффективность предложенного механизма распределения вычисле- 
ний в генетических алгоритмах. 


1 Распределенные вычисления в генетических алгоритмах 


Существует три основных стратегии использования распределенных вычислений 
в генетических алгоритмах [11]: 

— модель «мастер — раб» (тазет-Зауе СА$5). Особенность подхода состоит в син- 
хронной работе популяций. Популяции обладают общим адресным пространством. 
Мастер или управляющий процесс осуществляет развитие популяции, а подчинен- 
ные занимаются расчетом целевой функции; 

— модель «мелкоструктурного подчинения» (Ёпе-оташеа СА$5). Имеется управ- 
ляющий процесс и подчиненные сильно от него зависят, но работают в асинхронном 
режиме, при использовании общих ресурсов. В этой связи необходимо решать воп- 
росы синхронизации; 

— модель «сетевого взаимодействия» (соагзе-отате4 СА$5). Каждый процесс мо- 
делирует свою популяцию, используя собственное адресное пространство. Взаимо- 
действие с другими популяциями определяется из описания связей (топологии сети 
взаимодействия). 

Главное достоинство первой стратегии — это простота реализации, потому что 
распараллеливание происходит лишь на уровне вычисления ЦФ. Недостаток — низкая 
эффективность распараллеливания, большая нагрузка на управляющий процесс и 
сильная зависимость от него. 

Рассмотрим вторую стратегию на примере. В отличие от модели «мастер — раб» 
здесь появляется оператор миграции (обмен хромосомами между популяциями). 
Мелкоструктурное подчинение подразумевает миграцию хромосом только с соседни- 
ми популяциями (4 популяции для миграции). На первой итерации каждая популяция 
генерируется случайным образом, и ЦФ отличается от всех остальных популяций. 
В процессе работы ГА образуются области популяций со схожим значением ЦФ. По 
мере работы алгоритма эти области растут и конкурируют между собой. Такой меха- 
низм позволяет распространять «хорошие» индивиды на всей распределенной попу- 
ляции. Из недостатков второй стратегии можно выделить необходимость синхрони- 
зации, из-за асинхронного режима работы и общего адресного пространства. 

Модель сетевого взаимодействия часто сравнивают с сетью островов, которая 
рассматривает относительно изолированные общины, поэтому их называют «остров- 
ными» параллельными ГА [12]. Рассмотрим модель островов на примере. Пусть у 
нас есть 16 процессов и 1600 хромосом. Разобьем их на 16 отдельных популяций по 
100 особей. Каждая их них будет развиваться отдельно с помощью генетического 
алгоритма. Таким образом, можно сказать, что мы распределили особи по 16 изоли- 
рованным островам. Далее в зависимости от взаимосвязей островов осуществляется 
миграция. Оператор миграции можно выполнять как периодически, так и при на- 
ступлении определенного события (например, сходимости популяции к локальному 
оптимуму). Важно правильно подобрать частоту миграции (интервал времени между 
миграциями) и количество особей, участвующих в миграции. Частая миграция (как и 
миграция большого числа хромосом) приведет к смешению всех популяций и преж- 
девременной сходимости ГА. Редкая миграция также не сможет предотвратить преж- 
девременного схождения популяций. Островная модель позволяет каждому острову 
развиваться самостоятельно. Различные стратегии островов позволяют комбини- 
ровать различные архитектуры ГА и проводить исследования их взаимодействия. Из 
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недостатков можно выделить сложность реализации, из-за необходимости разработки 
механизмов взаимодействия отдельных островов (популяций). Стратегия сложна для 
понимания в силу того, что эффекты миграции не полностью понятны. Для расчета 
числа представителей схемы в одной из популяций необходимо определить модель 
оператора миграции. Если миграция будет происходить каждую итерацию (что слиш- 
ком часто, это приведет к вырождению популяций), то расчет можно вести по той же 
формуле, что и для модели мелкоструктурного подчинения. Если же миграция осу- 
ществляется по событию, то в формулу надо внести изменения, связанные с вероят- 
ностью наступления этого события: 


1-2 


где Р(М) -— это вероятность миграции, а К - число связных островов. Все это означает, 
что разработка моделей стратегии островов и их исследование — актуальная задача. 

В островной стратегии использования параллельных вычислений в генетических 
алгоритмах стратегии качество решения зависит от структуры связей, схемы миграции 
и многих других факторов. Островная модель мало изучена, особенно важное влияние 
на результаты имеет механизм миграции и подбор его параметров. В этой связи стра- 
тегия островов интересна для исследования и разработки различных модификаций. 


т(Н,1+ > т(Н,#)х Го» = року - ЦнухРОМ) (1) 


2 Нечеткий оператор миграции 


Работа М. Голубя и Д. Джакобовича подтверждает влияние частоты миграции 
на работу РГА [13]. Слишком частая миграция приводит к вырождению популяций и 
попаданию в локальный оптимум [14]. Для регулирования частоты миграции приме- 
няются различные методы, которые можно разбить на два типа: адаптивные и собы- 
тийные. В первом используются методы адаптации для настройки частоты миграции 
в процессе работы алгоритма. Во втором применяются событийные методы (мигра- 
ция осуществляется при наступлении определенных событий). Событием может быть 
уменьшение разнообразия в популяции ниже порогового значения или же низкая 
динамика поиска по сравнению с другими процессами. 


Рисунок 1 — Нечеткий контроллер для работы оператора миграции 
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Кратко опишем идею работы нечеткого оператора миграции: предлагается ис- 
пользовать показатель разнообразия в популяции и его сравнение с показателем раз- 
нообразия среди всех популяций. Для этого вводится нечеткое понятие «степень раз- 
нообразия популяции», представляемое нечетким множеством М = {низкая степень, 
высокая степень}. Далее формируется база правил, например, «если степень разно- 
образия низкая, то необходима миграция с 3 островами». На каждой итерации рас- 
считывается значение разнообразия, которое поступает в нечеткий логический конт- 
роллер, на выходе которого управляющее воздействие выполняет миграцию и если 
да, то с каким числом популяций (рис. 1). Преимущество метода заключается в воз- 
можности учета различных параметров популяции (кроме степени разнообразия). 
Отличие от других методов в том, что это не адаптация параметров работы алго- 
ритма, а реальная оценка эволюции в популяции. 

Степень разнообразия в популяции определяется как генотипическое расстоя- 
ние между всеми хромосомами. Если популяция стремится к вырождению (низкая 
степень разнообразия), то для предотвращения вырождения нужен оператор мигра- 
ции. Степень разнообразия будем рассчитывать по формуле: 


№ тп | | 
хх 


№ 
5 ра НА (2) 
4 2, М 


где №М — количество хромосом в популяции в ЦФ, и — число генов в хромосоме, 


х''-х/, — разница значения генов двух хромосом. Если хромосомы (альтернативные 
решения проблемы) в популяции отличаются друг от друга, значит, степень разнооб- 
разия в популяции достаточная для дальнейшей оптимизации. Кроме степени разно- 
образия в популяции надо учитывать уровень развития популяции, другими словами, 
отличие ЦФ в популяции от лучшего значения среди всех популяций. Обозначим 
этот критерий миграции как уровень развития и будем рассчитывать по формуле: 


Ави) 
тах(Ё,Р,,...Е‹)’ 
где /у — это ЦФ хромосомы в популяции, а Ёк — лучшее значение ЦФ среди всех 


популяций. Уровень развития представим нечетким множеством {низкий уровень 
развития, высокий уровень развития}. Составим базу правил (табл. 1). 


(3) 


4еу 


Таблица 1 — Нечеткие правила 


Уровень развития 
НИЗКИЙ высокий 

и низкая | Необходима | Необходима 
Е миграция, миграция, 
5 мигрирует | мигрирует 10% 
8 25% от от популяции 
Е популяции 
2. высокая | Необходима Нет 
ы миграция, | необходимости 
о 
= мигрирует в миграции 
= 0 
5 5% от 

популяции 
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3 Экспериментальные исследования 


Для исследования проведем эксперименты по оптимизации следующих функ- 
ций [16]: 
— максимизации функции Растригина [6]: 
Е(ж».. 20) = У (1008(2лх,-х2)-100}, (4) 


— минимизации функции Растригина (рис. 2): 
Е(х».. 0) = У. (х2 —1008(2лх,) +10), (5) 


в функциях Растригина переменная х; принимает значения в интервале [|-5.12; 5.12]; 
— минимизации функции сферы (рис. 3): 


Е(х, Г ЙЕ | )= 2. (6) 
— минимизации функции ее (рис. 4) 
м (7) 


Е(х,...ж)= 400 а. -Пе®| 1 


где х; принимает значения в интервале 600.0: 600.0]. 


-5 
5 


Рисунок 2 — Двумерная функция Рисунок 3 — Двумерная функция сферы 
Растригина 


Рисунок 4 — Двумерная функция Гриванкса 
Это сложные многоэкстремальные функции большой размерности (десять пе- 


ременных). Хромосома представляет собой вектор длины 10, содержащий значение 
переменных х....Хл. 
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Параметры ГА следующие: 

— размер популяции - 15; 

— число итераций -— 3000; 

— вероятность кроссинговера -— 75; 

— вероятность мутации - 5; 

— число процессов - 8. 

Для исследования эффективности работы нечеткого оператора миграции прове- 
дем его сравнение с простым оператором миграции. Под простым оператором мигра- 
ции понимается миграция через каждые п итераций по две хромосомы. Для наглядного 
исследования возьмем и равным 5, 10, 15, 20 и 25 итераций. Результат сравнения для 
задачи максимизации функции Растригина показан на графике (рис. 5). Для задачи 
минимизации функции Растригина - на рис. 6, и для задач минимизации функции 
сферы и функции Гриванкса на рис. 7 и 8 соответственно. 


Стратегии миграции 
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№ 25 итераций 


0 20 итераций 
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-10 -8 -6 4 -2 0 
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Рисунок 5 — График сравнения результатов работы алгоритма с различными 
моделями оператора миграции для задачи максимизации функции Растригина 
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Рисунок 6 — График сравнения результатов работы алгоритма с различными 
моделями оператора миграции для задачи минимизации функции Растригина 
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Рисунок 7 — График сравнения результатов работы алгоритма с различными 
моделями оператора миграции для задачи минимизации функции сферы 
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Рисунок 8 — График сравнения результатов работы алгоритма с различными 
моделями оператора миграции для задачи минимизации функции Гриванкса 


Все исследования проводились в асинхронном режиме, на компьютере с дву- 
ядерным процессором фирмы АМЮ, с объемом оперативной памяти в 1Гб, под управ- 
лением операционной системы \/ш@до\з ХР. Результаты исследований показали 
возможность эффективного использования РГА на многоядерных процессорах. 
Исследования показали, что с увеличением размера буфера хромосом решение улуч- 
шается, но наблюдается так называемая точка насыщения, когда дальнейшее увели- 
чение объема буфера будет не целесообразно. В этой связи предложенная модель по- 
зволяет проводить исследования механизма миграции как одной из основных состав- 
ляющих распределенных вычислений в ГА. Исследования показали, что необходимо 
использовать динамические связи между популяциями. Гибридный оператор мигра- 
ции позволяет проводить миграцию только когда это необходимо, например, при 
угрозе преждевременной сходимости. 


Выводы 


В настоящий момент доступны крупные вычислительные центры, так называемые 
мейнфреймы. Все они работают в многопоточном режиме. В этой связи разработка 
моделей использования параллельных вычислений в многопоточной среде — актуаль- 
ная задача, т.к. развитие многопроцессорных вычислительных машин продолжается, и 
для их эффективного использования необходимы алгоритмы распределенных вычи- 
слений. 

Предложенная модель позволила отойти от жестко заданных связей между про- 
цессами и показала эффективность динамических связей между ними. Модифициро- 
ванная модель доказала важность оператора миграции и влияние механизма опера- 
тора миграции на работу РГА. Для развития РГА необходима разработка нового 
механизма миграции, использующего динамические связи процессов. 
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Д.С. Книш, В.М. Курейчик 

Паралельний генетичний алгоритм з нечиким оператором миграций 

Розвиток обчислювально! технйки сприяв розвитку метод1в паралельних обчислень при розв’язанн! 
оптимзащйних задач. Генетичн! алгоритми — ефективний зас1б розв’язання р1зних оптимзащйних 
задач. У зв’язку з цим необхдн! стратегИ застосування паралельних обчислень в генетичних алгоритмах. 
За аналогею з еволющйними процесами за основу паралельних обчислень у генетичних алгоритмах 
взято механ1зм миграций. Авторами пропонуеться реалзащя оператора миграцй з використанням апарату 
нечтко! логки. Проведен! випробування розроблено1 модел1 показали Й ефективнисть. 


р.5. Кпузй, Г.М. КугасШЕ 

Тве Обида Сепейс Аюогт ув Еи77у МогаНоп Орегаог 

ТБе сотрийпе епошеетте 4еуе!орше Ваз а$$15{е4 песеззНу оЁ гезеагсВ аррИсаНопз о# 415616 \е4 сотринпе 
тео@$ ЕЮюг орнпитайоп 1азК$. ТВе сепейс а1оогит 15 еЙесйуе шзилитет Гог Ч1егеп! орйпитайоп {а5К$ 
зо[уш?. ТВегеироп те фо@$ оЁ аррПсайоп @зитощед сотрийпе аге гедитеа Бу сепейс а]еогИллз. ТВе 
рагаПе| сотрийп? ш сепейс а1еоги лаз Базе оп пиотайоп орегаюг Бу апа!огу Ул еуоаопагу ргосез$ 
ш пабхге. ТБе тесвапзт оЁ пизтайоп \ИБ Ё77у 1021с изшо 1$ зассезе4. Тез оЁ пиотайоп орегаюг Ваз 
зво\ме4 и еЙесйуепезз. 


Статья поступила в редакцию 15.07.2010. 
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